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Using 2.93 fb! of data taken at 3.773 GeV with the BESIII detector operated at the BEPCII collider, we 
study the semileptonic decays D* — K?e* v, and Dt — n?etve. We measure the absolute decay branching 
fractions B(D* > K?e*v.) = (8.60 + 0.06 + 0.15) x 10? and B(D* — ret.) = (3.63 + 0.08 + 
0.05) x 107°, where the first uncertainties are statistical and the second systematic. We also measure the 
differential decay rates and study the form factors of these two decays. With the values of |V..| and |V-a| from 
Particle Data Group fits assuming CKM unitarity, we obtain the values of the form factors at q^ = 0, fi (0) = 
0.725 +0.004+ 0.012 and f£ (0) = 0.622 + 0.012 + 0.003. Taking input from recent lattice QCD calculations 
of these form factors, we determine values of the CKM matrix elements |V-s| = 0.944 + 0.005 + 0.015 + 0.024 


and |Vea| = 0.210 + 0.004 + 0.001 + 0.009, where the third uncertainties are theoretical. 
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PACS numbers: 13.20.Fc, 12.15.Hh 


I. INTRODUCTION 


In the Standard Model (SM) of particle physics, the mixing 
between the quark flavours in the weak interaction is param- 
eterized by the Cabibbo-Kobayashi-Maskawa (CKM) matrix, 
which is a 3 x 3 unitary matrix. Since the CKM matrix el- 
ements are fundamental parameters of the SM, precise deter- 
minations of these elements are very important for tests of the 
SM and searches for New Physics (NP) beyond the SM. 

Since the effects of strong and weak interactions can be well 
separated in semileptonic D decays, these decays are excel- 
lent processes from which we can determine the magnitude 
of the CKM matrix element V..(q). In the SM, neglecting the 
lepton mass, the differential decay rate for Dt — Petr, 
(P = K? or 7°) is given by [1] 


dU G2 

E = Xj Vew Pf (P, (1) 
where G p is the Fermi constant, Ves(a) is the corresponding 
CKM matrix element, p is the momentum of the meson P in 
the rest frame of the D meson, q? is the squared four momen- 
tum transfer, i.e., the invariant mass of the lepton and neutrino 
system, and f, (q?) is the form factor which parameterizes 
the effect of the strong interaction. In Eq. (1), X is a multi- 
plicative factor due to isospin, which equals to 1 for the decay 
Dt — K°ety, and 1/2 for the decay Dt — n? e*v,. 

In this article, we report the experimental study of D^ — 
K?etv, and Dt — a? e*v, decays using a 2.93 fb! [2] data 
set collected at a center-of-mass energy of \/s = 3.773 GeV 
with the BESIII detector operated at the BEPCII collider. 
Throughout this paper, the inclusion of charge conjugate 
channels is implied. 

The paper is structured as follows. We briefly describe 
the BESIII detector and the Monte Carlo (MC) simulation 
in Sec. II. The event selection is presented in Sec. III. The 
measurements of the absolute branching fractions and the dif- 
ferential decay rates are described in Sec. IV and V, respec- 
tively. In Sec. VI we discuss the determination of form factors 
from the measurements of decay rates, and finally, in Sec. VII, 
we present the determination of the magnitudes of the CKM 
matrix elements V., and V.g. A brief summary is given in 
Sec. VIII. 


Il. BESIII DETECTOR 


The BESIII detector is a cylindrical detector with a solid- 
angle coverage of 93% of 47, designed for the study of hadron 
spectroscopy and 7-charm physics. The BESIII detector is de- 
scribed in detail in Ref. [3]. Detector components particularly 
relevant for this work are (1) the main drift chamber (MDC) 
with 43 layers surrounding the beam pipe, which performs 
precise determination of charged particle trajectories and pro- 
vides a measurement of the specific ionization energy loss 
(dE/dx); (2) a time-of-flight system (TOF) made of plastic 


scintillator counters, which are located outside of the MDC 
and provide additional charged particle identification informa- 
tion; and (3) the electromagnetic calorimeter (EMC) consist- 
ing of 6240 CsI(TI) crystals, used to measure the energy of 
photons and to identify electrons. 

A GEANT4-based [4] MC simulation software [5], which 
contains the detector geometry description and the detector re- 
sponse, is used to optimize the event selection criteria, study 
possible backgrounds, and determine the reconstruction effi- 
ciencies. The production of the (3770), initial state radiation 
production of 1/(3686) and J/1/, as well as the continuum pro- 
cesses of e* e^ — r*vr- ande*e^ — qq (q = u,d,s) are 
simulated by the MC event generator KKMC [6]; the known 
decay modes are generated by EVTGEN [7] with the branching 
fractions set to the world average values from the Particle Data 
Group (PDG) [8]; while the remaining unknown decay modes 
are modeled by LUNDCHARM [9]. We also generate signal 
MC events consisting of (3770) — D+ D- events in which 
the D~ meson decays to all possible final states and the D+ 
meson decays to a hadronic or a semileptonic decay final state 
being investigated. In the generation of signal MC events, the 
semileptonic decays D^ — K?e*v, and Dt — «Petty, 
are modeled by the the modified pole parametrization (see 
Sec. VI A). 


III. EVENT RECONSTRUCTION 


The center-of-mass energy of 3.773 GeV corresponds to the 
peak of the (3770) resonance, which decays predominantly 
into DD (D? D? or D* D-) meson pairs. In events where 
a D' meson is fully reconstructed, the remaining particles 
must all be decay products of the accompanying D^ meson. 
In the following, the reconstructed meson is called “tagged 
D=” or “D~ tag". In a tagged D~ data sample, the recoiling 
D* decays to K°etv, or neve can be cleanly isolated and 
used to measure the branching fraction and differential decay 
rates. 


A. Selection of D~ tags 


We reconstruct D~ tags in the following nine hadronic 
modes: D + K*» * ,D > Klr, DT > KEK-, 
D- > Kt K-1-, D7 K'a nr nrl, D ^ ata |, 
D^ > K8n-7, D- — K*a-« "n rnt, and D- > 
Kea aa. The selection criteria of D~ tags used here 
are the same as those described in Ref. [10]. Tagged D^ 
mesons are identified by their beam-energy-constrained mass 


Mpo =, [ER 6 — |Piag|?/c?, where Eyecam is the beam 


1 We veto the K oa candidates when a + «- invariant mass falls within 
the K 9 mass window. 
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energy, and tag is the measured 3-momentum of the tag can- 
didate ?. We also use the variable AE — Erag — P'beam; 
where Frag is the measured energy of the tag candidate, to 
select the D~ tags. Each tag candidate is subjected to a tag 
mode-dependent A E requirement as shown in Table I. If there 
are multiple candidates per tag mode for an event, the one with 
the smallest value of |AF| is retained. 

The Msc distributions for the nine D~ tag modes are 
shown in Fig. 1. A binned extended maximum likelihood fit 
is used to determine the number of tagged D~ events for each 
of the nine modes. We use the MC simulated signal shape 
convolved with a double-Gaussian resolution function to rep- 
resent the beam-energy-constrained mass signal for the D7 
daughter particles, and an ARGUS function [11] multiplied by 
a third-order polynomial [12, 13] to describe the background 
shape for the Mgc distributions. In the fits all parameters of 
the double-Gaussian function, the ARGUS function, and the 
polynomial function are left free. The solid lines in Fig. 1 
show the best fits, while the dashed lines show the fitted back- 
ground shapes. The numbers of the D~ tags (Vtag) within the 
Mgc signal regions given by the two vertical lines in Fig. 1 
are summarized in Table I. In total, we find 1703054 + 3405 
single D~ tags reconstructed in data. The reconstruction ef- 
ficiencies of the single D~ tags, €¢ag, as determined with the 
MC simulation, are shown in Table I. 


B. Reconstruction of semileptonic decays 


Candidates for semileptonic decays are selected from the 
remaining tracks in the system recoiling against the D~ tags. 
The dE/dx, TOF and EMC measurements (deposited energy 
and shape of the electromagnetic shower) are combined to 
form confidence levels for the e hypothesis (C L.), the m hy- 
pothesis (CL), and the K hypothesis (C Li). Positron can- 
didates are required to have CLe greater than 0.1% and to 
satisfy CLe/(CLe + CLr +CLx) > 0.8. In addition, we 
include the 4-momenta of near-by photons within 5? of the 
direction of the positron momentum to partially account for 
final-state-radiation energy losses (FSR recovery). The neu- 
tral kaon candidates are built from pairs of oppositely charged 
tracks that are assumed to be pions. For each pair of charged 
tracks, a vertex fit is performed and the resulting track param- 
eters are used to calculate the invariant mass, M (z*«- ). If 
M (1*7) is in the range (0.484, 0.512) GeV/c?, the mt «^ 
pair is treated as a K$ candidate and is used for further anal- 
ysis. The neutral pion candidates are reconstructed via the 
n? — yy decays. For the photon selection, we require the 
energy of the shower deposited in the barrel (end-cap) EMC 
greater than 25 (50) MeV and the shower time be within 
700 ns of the event start time. In addition, the angle be- 
tween the photon and the nearest charged track is required 
to be greater than 10°. We accept the pair of photons as a 7° 


? Tn this analysis, all four-momentum vectors measured in the laboratory 
frame are boosted to the e^ e^ center-of-mass frame. 


candidate if the invariant mass of the two photons, M (^), 
is in the range (0.110, 0.150) GeV /c?. A 1-Constraint (1-C) 
kinematic fit is then performed to constrain M (yy) to the 7° 
nominal mass, and the resulting 4-momentum of the candidate 
n? is used for further analysis. 

We reconstruct the D+ — K?etv, decay by requiring ex- 
actly three additional charged tracks in the rest of the event. 
One track with charge opposite to that of the D~ tag is iden- 
tified as a positron using the criteria mentioned above, while 
the other two oppositely charged tracks form a Iv T candidate. 
For the selection of the D+ — meve decay, we require that 
there is only one additional charged track consistent with the 
positron identification criteria and at least two photons that are 
used to form a 7° candidate in the rest of the event. If there are 
multiple 7° candidates, the one with the minimum x? from the 
1-C kinematic fit is retained. In order to additionally suppress 
background due to wrongly reconstructed or background pho- 
tons, the semileptonic candidate is further required to have the 
maximum energy of any of the unused photons, Fy max, less 
than 300 MeV. 

Since the neutrino is undetected, the kinematic variable 
Umiss = Emiss — C|Pmiss| is used to obtain the information 
about the missing neutrino, where Emiss and piss are, respec- 
tively, the total missing energy and momentum in the event. 
The missing energy is computed from Emiss = Epeam — EP — 
E+, where Ep and E,. are the measured energies of the 
pseudoscalar meson and the positron, respectively. The miss- 
ing momentum Pmiss is given by Dmiss = Pp+ = PP m Ders 
where pp+, pp and p,+ are the 3-momenta of the D* me- 
son, the pseudoscalar meson and the positron, respectively. 
The 3-momentum of the D* meson is taken as pp+ = 
—Ptag y (Ebeam/C)? — (Mp+c)?, where fig is the direction 
of the momentum of the single D~ tag, and mp+ is the D* 
mass. If the daughter particles from a semileptonic decay are 
correctly identified, Umiss is near zero, since only one neutrino 
is missing. 

Figure 2 shows the Umiss distributions for the semileptonic 
candidates, where the potential backgrounds arise from the 
D D processes other than signal, (3770) — non-D D decays, 
ete — T*r-, continuum light hadron production, initial 
state radiation return to J/w and (3686). The background 
for D+ — K°etv, is dominated by D+ + K*(892)9e*v. 
and D+ — K? utv. For D+ — s? e*v,, the background is 
mainly from Dt > K}? etve and Dt > KQ(x?n?)e*v.. 

Following the same procedure described in Ref. [13], we 
perform a binned extended maximum likelihood fit to the 
Umiss distribution for each channel to separate the signal 
from the background component. The signal shape is con- 
structed from a convolution of a MC determined distribution 
and a Gaussian function that accounts for the difference of 
the Umiss resolutions between data and MC simulation. The 
background shape is formed from MC simulation. From the 
fits shown as the overlaid curves in Fig. 2, we obtain the yields 
of the observed signal events to be Nops(D+ > K?etve) = 
26008 + 168 and Nops(Dt > s?e* v.) = 3402 + 70, respec- 
tively. 

To check the quality of the MC simulation, we examine the 
distributions of the reconstructed kinematic variables. Fig- 


TABLE I. The AE requirements, the Mpc signal regions, the yields of the D' tags (Ntag) reconstructed in data, and the reconstruction 
efficiency (€tag) of D' tags. The uncertainties are statistical only. 


Tag mode AE (MeV) Mpc (GeV/c?) Nag Etag (96) 
D-— Kta-7 —45, 45) 1.8640, 1.8770) 806830 + 1070 51.8 £0.1 
D- > Ken7 —45, 45) 1.8640, 1.8770) 102755 + 372 56.2 + 0.2 
D- > KK- —45, 45) 1.8650, 1.8770) 19566 + 185 52.1 +0.5 
D > KK m —50, 50) 1.8650, 1.8780) 68216 + 966 41.2 +0.3 
D > Kronan’ —78, 78) 1.8620, 1.8790) 271571 + 2367 27.3 +0.1 
D > ntn r —45, 45) 1.8640, 1.8770) 32150 + 371 56.9 + 0.7 
D- > Ken 7° —T5, 75) 1.8640, 1.8790) 245303 + 1273 31.3 +0.1 
D > Ktn on nr n —52, 52) 1.8630, 1.8775) 30923 + 733 22.1 + 0.2 
D + Kgn n rt —50, 50) 1.8640, 1.8770) 125740 + 1203 33.0 + 0.2 
Sum 1703054 + 3405 
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FIG. 1. Fits (solid lines) to the Mpc distributions (points with error bars) in data for nine D' tag modes. The two vertical lines show the 
tagged D' mass regions. 


~ x10? 3 x10? ure 3 shows the comparisons of the momentum distributions 


(b) | of data and MC simulation. 


IV. BRANCHING FRACTION MEASUREMENTS 


A. Determinations of branching fractions 


-0.2 0 0.2 0.4 


The branching fraction of the semileptonic decay D* — 


Uriss (GeV) Pe*v, is obtained from 
= + + 
FIG. 2. Distributions of Umiss for the selected (a) Dt + K?’ etve B(D* > Petve) = Nass (D* — Pe*ve) , 2) 
and (b) D? — a? e*v, candidates (points with error bars) with fit i Niage(Dt > Petve) 
projections overlaid (solid lines). The dashed curves show the back- . 
ground determined by the fit. where N;4, is the number of D~ tags (see Sec. IMA), 


Nops(D* — Petve) is the number of observed Dt — 
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FIG. 3. Momentum distributions of selected events (with |Umiss| < 
60 MeV) for (a) K?, (b)e* from D* — K°et ve, (c) 7°, and (d) e* 
from D+ — a? e*v,. The points with error bars represent data, the 
(blue) open histograms are MC simulated signal plus background, 
the shaded histograms are MC simulated background only. 


Pe*vw, decays within the D~ tags (see Sec. IIIB), and 
e(D* — Petv,) is the reconstruction efficiency. Here the 
D+ — Ke*v, efficiency includes the K fraction of the 
K? and KÌ — ntr” branching fraction, the Dt > net v, 
efficiency includes the 7° — yy branching fraction [8]. 

Due to the difference in the multiplicity, the reconstruc- 
tion efficiency varies slightly with the tag mode. For each 
tag mode i, the reconstruction efficiency is given by e^ = 
Etag SL / tabs where the efficiency for simultaneously find- 
ing the D — Petve semileptonic decay and the D^ me- 
son tagged with mode 7, Clog SLs is determined using the sig- 
nal MC sample, and Shag is the corresponding tag efficiency 
shown in Table I. These efficiencies are listed in Table II. The 
reconstruction efficiency for each tag mode is then weighted 
according to the corresponding tag yield in data to obtain the 
average reconstruction efficiency, € = }0;(Nfagé")/Ntag, as 
listed in the last row in Table II. 

Using the control samples selected from Bhabha scatter- 
ing and DD events, we find that there are small discrepan- 
cies between data and MC simulation in the positron track- 
ing efficiency, positron identification efficiency, K9 and 7° 
reconstruction efficiencies. We correct for these differences 
by multiplying the raw efficiencies e(D+ — K?etve) and 
e(Dt — netve) determined in MC simulation by factors 
of 0.9957 and 0.9910, respectively. The corrected efficiencies 
are found to be e''(D* — K?etve) = (17.75 + 0.03)% and 
Ee (Dt — n°et vy.) = (55.02 + 0.10)%, where the uncertain- 
ties are only statistical. 

Inserting the corresponding numbers into Eq. (2) yields the 
absolute decay branching fractions 


B(D* — K?e*ve) = (8.60 + 0.06 


0.15) x 107? (3) 


and 


B(D* > a?9e*v,) = (3.63 + 0.08 


0.05) x 107%, (4 


where the first uncertainties are statistical and the second sys- 
tematic. 


B. Systematic uncertainties 


The systematic uncertainties in the measured branching 
fractions of D+ — K?e*v, and Dt — nPe*v, decays in- 
clude the following contributions. 

Number of D~ tags. The systematic uncertainty of the 
number of D~ tags is 0.5% [10]. 

e* tracking efficiency. Using the positron samples selected 
from radiative Bhabha scattering events, the e* tracking effi- 
ciencies are measured in data and MC simulation. Consider- 
ing both the polar angle and momentum distributions of the 
positrons in the semileptonc decays, a correction factor of 
1.0021 + 0.0019 (1.0011 + 0.0015) is determined for the e+ 
tracking efficiency in the branching fraction measurement of 
D* + K?etw, (Dt — n?etv,) decay. This correction is 
applied and an uncertainty of 0.19% (0.15%) is used as the 
corresponding systematic uncertainty. 

e* identification efficiency. Using the positron samples se- 
lected from radiative Bhabha scattering events, we measure 
the e* identification efficiencies in data and MC simulation. 
Taking both the polar angle and momentum distributions of 
the positrons in the semileptonic decays into account, a cor- 
rection factor of 0.9993 + 0.0016 (0.9984 + 0.0014) is deter- 
mined for the e* identification efficiency in the measurement 
of B(D* — K?etve) (B(D* — nPe*v,)). This correction 
is applied, and an amount of 0.16% (0.14%) is assigned as the 
corresponding systematic uncertainty. 

KÌ and m? reconstruction efficiency. The momentum- 
dependent efficiencies for K9 (19) reconstruction in data and 
in MC simulation are measured with DD events. Weighting 
these efficiencies according to the K E: (1?) momentum dis- 
tribution in the semileptonic decay leads to a difference of 
(—0.57 + 1.62)96 ((—0.85 + 1.00)%) between the K9 (x?) 
reconstruction efficiencies in data and MC simulation. Since 
we correct for the systematic shift, the uncertainty of the cor- 
rection factor, 1.62% (1.00%), is taken as the corresponding 
systematic uncertainty in the measured branching fraction of 
Dt — K°ety, (Dt — wetv.). 

Requirement on Ey max. By comparing doubly tagged DD 
hadronic decay events in the data and MC simulation, the sys- 
tematic uncertainty due to this source is estimated to be 0.1%. 

Fit to the Umiss distribution. To estimate the uncertain- 
ties due to the fits to the Umiss distributions, we refit the 
Umiss distributions by varying the bin size and the tail pa- 
rameters (which are used to describe the signal shapes and 
are determined from MC simulation) to obtain the number of 
signal events from Dt semileptonic decays. We then com- 
bine the changes in the yields in quadrature to obtain the sys- 
tematic uncertainty (0.12% for D+ > K?etv,, 0.52% for 
D* — n?etve). Since the background function is formed 
from many background modes with fixed relative normaliza- 
tions, we also vary the relative contributions of several of the 
largest background modes based on the uncertainties in their 
branching fractions (0.12% for D+ — K°etv,, 0.01% for 


TABLEII. The reconstruction efficiencies for D+ > K°etv. and Dt — n?e” v, determined from MC simulation. The efficiencies include 
the branching fractions for K° and 7°. The uncertainties are statistical only. 


Tag mode Etag,sL(Dt — K'e*wv,) (9) e(D* — K?e*v.) (9) eugsL(D* — mete) (9) e(D* — n? e* v.) (96) 
D —K'nm« 9.21 + 0.02 17.77 £0.04 28.44 + 0.06 54.88 £0.13 
D- > K3n7 10.14 + 0.05 18.05 + 0.11 31.15 £0.15 55.43 + 0.34 
D- > KOK- 9.30 + 0.08 17.84 + 0.22 28.68 + 0.23 55.02 + 0.67 
D- > Kt K-17 7.39 + 0.06 17.92 + 0.18 22.53 + 0.16 54.66 + 0.53 
D-— Kn «x? 4.98 + 0.02 18.25 + 0.09 15.49 + 0.06 56.72 + 0.29 
D^ > ntron 10.44 + 0.11 18.34 0.30 32.93 + 0.33 57.82 + 0.94 
D~ + Kon m? 5.67 + 0.01 18.11 + 0.08 17.83 + 0.04 56.92 + 0.25 
D-— Kn mn n n? 3.50 + 0.04 15.88 + 0.25 11.74 + 0.14 53.20 £0.81 
D > Kin zv nt 5.55 + 0.02 16.84 + 0.14 18.12 + 0.06 54.97 + 0.45 
Average 17.83 + 0.03 


55.52 + 0.10 


D* — mÜe*y,). In addition, we convolute the background 
shapes formed from MC simulation with the same Gaussian 
function in the fits (0.02% for Dt — K°etv,, 0.30% for 
D+ — n?e*v,). Finally we assign the relative uncertainties 
to be 0.2% and 0.696 for Dt + K?e*v, and Dt > m?etve, 


TABLE III. Summary of the systematic uncertainties considered in 
the measurements of the branching fractions of D+ + Ke*v, and 
Dt + metre decays. 


Systematic uncertainty (%) 


respectively. Source D H — K°e Ve D H — re H Ve 
Form factor. In order to estimate the systematic uncertainty Number of D~ tags 0.5 0.5 
: : : Tracking for e* 0.19 0.15 
associated with the form factor used to generate signal events Pn 
in the MC simulation, we re-weight the signal MC events so PID EE ; 0:16 0.14 
2 . K's reconstruction 1.62 sai 
that the q^ spectra agree with the measured spectra. We then o : 
: i : um 7T reconstruction ens 1.00 
remeasure the branching fraction (partial decay rates in dif- Reaui 
uas : , ] . equirement on Fy max 0.1 0.1 
ferent q^ bins) with the newly weighted efficiency (efficiency Fit to Umiss distribution 0.2 0.6 
matrix). The maximum relative change of the branching frac- Form factor 0.2 0.2 
tion (partial decay rates in different q? bins) is 0.2% and is FSR recovery 0.1 0.5 
assigned as the systematic uncertainty. MC statistics 0.2 0.2 
FSR recovery. The differences between the results with KÌ /n? branching fraction 0.07 0.03 
FSR recovery and the ones without FSR recovery are assigned Total 1.76 1.41 


as the systematic uncertainties due to FSR recovery. We find 
the differences are 0.1% and 0.5% for Dt — K?etve and 
Dt — ret, respectively. 

MC statistics. The uncertainties in the measured branching 
fractions due to the MC statistics are the statistical fluctua- 
tion of the MC samples, which are 0.2% for both of D* > 
K?ety, and D+ — 7°e+v, semileptonic decays. 

K} and n? decay branching fractions. We include an un- 
certainty of 0.0796 (0.0396) on the branching fraction mea- 
surement of Dt > K?e*v, (Dt — met.) to account 
for the uncertainty of the branching fraction of K} — ntr- 
(19 > yy) decay [8]. 

Table III summarizes the systematic uncertainties in the 
measurement of the branching fractions. Adding all system- 


shown in Fig. 4. Our measured branching fractions are in 
agreement with the other experimental measurements, but are 
more precise. For Dt — wet v,, our result is lower than the 
only other existing measurement by CLEO-c [15] by 2.0c. 


Using our previous measurements of B(D? + K~etv-) 
and B(D? — «- etve) [13], the results obtained in this anal- 
ysis, and the lifetimes of D? and D+ mesons [8], we obtain 
the ratios 


0 x ulli 
atic uncertainties in quadrature yields the total systematic un- Igz Maus AO =1.0340.0140.02 (5) 
certainties of 1.76% and 1.41% for Dt — K?etve and [( Dt  K'etu) 
Dt — n? e*v,, respectively. 
and 
C. Comparison 
T(D9? > metre 
pa T e^") 1034 0.03 0.02, (6) 


The comparisons of our measured branching fractions for 
D+ + K?etv and D+ — a? e*v, decays with those pre- 
viously measured at the BES-II [14], CLEO-c [15] and BE- 
SIII [16, 17] experiments as well as the PDG values [8] are 


2T(D* — n? etu.) 


which are consistent with isospin symmetry. 
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BES-II[14]  8.95+1.59+0.67 

CLEO-c [15] 8.83+0.10+0.20 

PDG [8] 8.90£0.15 

BESIII [16]  8.962+0.054+0.206 

(KK) 

BESIII [17]  8.59+0.14+0.21 
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FIG. 4. 
correspond to the 1c limits of the world averages. 


V. PARTIAL DECAY RATE MEASUREMENTS 
A. Determinations of partial decay rates 


To study the differential decay rates, we divide the semilep- 
tonic candidates satisfying the selection criteria described in 
Sec. III into bins of g?. Nine (seven) bins are used for 
Dt + K?etw, (Dt — m?etv,). The range of each bin 
is given in Table IV. The squared four momentum transfer 
q? is determined for each semileptonic candidate by q? = 
(Ee + Ey, )? [ c* — (B. +B)? / ?, where the energy and mo- 
mentum of the missing neutrino are taken to be E, = Emiss 
and py, =  Ewisfwiss/C, respectively. For each q^ bin, 
we perform a maximum likelihood fit to the corresponding 
Umiss distribution following the same procedure described in 
Sec. IIL B and obtain the signal yields as shown in Table IV. 

To account for detection efficiency and detector resolution, 
the number of events N?,. observed in the ith q? bin is ex- 
tracted from the relation 


Nbins 

A y j j 

obs ^ ej Nora? (7) 
j-1 


where Npins is the number of q? bins, NÍ, q is the number of 


semileptonic decay events produced in the tagged D~ sample 
with the q? filled in the jth bin, and e; j is the overall efficiency 
matrix that describes the efficiency and smearing across q? 
bins. The efficiency matrix element &;; is obtained by 


nie 1 

J 

Eij = —gm —— fij, (8) 
nj Etag 


where n;** is the number of the signal MC events generated 


ij 
in the jth q? bin and reconstructed in the ith q? bin, n°" is 
the total number of the signal MC events which are generated 
in the jth q? bin, and fij is the matrix to correct for data- 
MC differences in the efficiencies for e* tracking, e* iden- 


tification, and K? (7°) reconstruction. Table V presents the 


Comparison of the branching fraction measurements for D^ — K°etv, (left) and Dt > m 


CLEO-c [15] 4.05+0.16+0.09 

PDG [8] 4.05+0.18 

This work 3.63+0.08+0.05 —— 
phe te Eraa aaa laa baaalaanla 
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?e*y, (right). The green bands 


average overall efficiency matrices for D+ > K?etv, and 
D* — m«Pe*v, decays. To produce this average overall ef- 
ficiency matrix, we combine the efficiency matrices for each 
tag mode weighted by its yield shown in Table I. The diag- 
onal elements of the matrix give the overall efficiencies for 
D* — Pe*v, decays to be reconstructed in the correct q? 
bins in the recoil of the single D~ tags, while the neighboring 
off-diagonal elements of the matrix give the overall efficien- 
cies for cross feed between different q? bins. 

The partial decay width in the ith bin is obtained by invert- 
ing the matrix Eq. (7), 


AT; = Nha - 1 
‘ Tp+ Mag 


Tp Ni tag 
where 7p. is the lifetime of the D* meson [8]. The q?- 
dependent partial widths for Dt — K?e*v, and Dt > 
TÜe* y, are summarized in Table VI. Also shown in Table VI 
are the statistical uncertainties and the associated correlation 
matrices. 


B. Systematic covariance matrices 


For each source of systematic uncertainty in the measure- 
ments of partial decay rates, we construct an Npins X Nbins 
systematic covariance matrix. A brief description of each con- 
tribution follows. 

D* lifetime. The systematic uncertainty associated with the 
lifetime of the D^ meson (0.7%) [8] is fully correlated across 
q? bins. 

Number of D~ tags. The systematic uncertainty from the 
number of the single D~ tags (0.5%) is fully correlated be- 
tween q? bins. 

et, K d and n? reconstruction. The covariance matri- 
ces for the systematic uncertainties associated with the e^ 
tracking, e* identification, AK, and 7° reconstruction effi- 
ciencies are obtained in the following way. We first vary 


TABLE IV. Summary of the range of each q? bin, the number of the observed signal events for Dt + K'e* v, and Dt — mlet v, in data. 


Dt + K°ety. 


Bin No. 1 2 3 4 5 6 7 8 9 

q? (GeV? /c*)  [0.0,0.2)  (0.2,0.4)  ([0.4,0.6) —([0.6,0.8)  [0.8,1.0)  [L0,1.2)  [1.2,1.4)  [L4,1.6)  [1.6, grax) 

Nobs 5842+81 4935 +73 4180+67 3515+62 2818+55 2120+48 1460+40 8604-31 302+19 
Dt neti 

Bin No. 1 2 3 4 5 6 7 

q?(GeV?/c*)  [0.0,0.3)  [0.3,0.6)  [0.6,0.9)  (0.9,1.2)  [12,15)  [1.5,2.0)  [2.0, q2ax) 

Nobs 6582-29 562+27 4672-25 448424 401+24 470+26 4042330 


TABLE V. Efficiency matrices £i; given in percent for Dt + K?e*v, and Dt > m?e*v, decays. The column gives the true q° bin j, 
while the row gives the reconstructed q? bin i. The statistical uncertainties in the least significant digits are given in the parentheses. 


D* + K? etve 


Rec. q True q? (GeV? /c’) 

(GeV?/c*)  [0.0,0.2)  [0.2,0.4)  [0.4,0.6)  ([0.6,0.8 [08,1.0  [1.0,1.2)  [1.2,1.4)  [L4,16)  [1L.6,924.) 

0.0, 0.2) 18.53(6) 0.95(1) 0.07(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 

0.2, 0.4) 0.37(1) 16.86(6) 1.03(2) 0.05(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 

0.4, 0.6) 0.00(0) 0.40(1) 16.03(6) 1.03(2) 0.03(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 

0.6, 0.8) 0.00(0) 0.00(0) 0.46(1) 15.72(6) 0.95(2) 0.02(0) 0.00(0) 0.00(0) 0.00(0) 

0.8, 1.0) 0.00(0) 0.00(0) 0.01(0) 0.44(1) 15.78(7) 0.93(2) 0.01(0) 0.00(0) 0.00(0) 

1.0, 1.2) 0.00(0) 0.00(0) 0.00(0) 0.01(0) 0.46(1) 15.76(8) 0.80(2) 0.01(0) 0.00(0) 

1.2, 1.4) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.42(1) 15.58(9) 0.74(3) 0.00(0) 

1.4, 1.6) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.38(2) 15.45(12)  0.78(5) 

1.6,424.) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.00(0) 0.28(2) 15.98(19) 
Dt + metre 

Rec. q? True q? (GeV? /c’) 

(GeV?/c?)  [0.0,0.3) | [0.3,0.6) | ([0.6,0.9 ^ [0.9,1.2)  [1.2,1.5)  [1.5,20) [2.0, q2ax) 

0.0, 0.3) 53.84(15) — 2.27(3) 0.17(1) 0.01(0) 0.00(0) 0.00(0) 0.00(0) 

0.3, 0.6) 4.00(5)  48.24(15)  2.31(4) 0.14(1) 0.00(0) 0.00(0) 0.01(0) 

0.6, 0.9) 0.14(1) 5.66(6) ^ 46.15(15)  2.34(4) 0.10(1) 0.00(0) 0.00(0) 

0.9, 1.2) 0.04(0) 0.22(1) 6.24(6) ^ 44.51(16) — 2.16(4) 0.05(0) 0.00(0) 

1.2, 1.5) 0.04(0) 0.08(1) 0.31(1) 6.33(7) | 43.33(17) 1.36(3) 0.02(0) 

1.5, 2.0) 0.03(0) 0.08(1) 0.22(1) 0.58(2) 6.52(8) ^ 45.48(16) 1.12(3) 

2.0, ghax)  0.13(1) 0.21 (1) 0.34(1) 0.68(2) 1.30(3) 5.52(6)  50.46(19) 


the corresponding correction factors according to their uncer- 
tainties, then remeasure the partial decay rates using the effi- 
ciency matrices determined from the re-corrected signal MC 
events. The covariance matrix due to this source is assigned 
via Ci; = Ó(AT;)Ó(AT ;), where 9(AT';) denotes the change 
in the partial decay rate measurement in the ith q? bin. 

Requirement on Ey max. We take the systematic uncer- 
tainty of 0.1% due to the E~ max requirement on the selected 
events in each q? bin, and assume that this uncertainty is fully 
correlated between q? bins. 

Fit to the Umiss distribution. The technique of fitting the 
Umiss distributions affects the number of signal events ob- 
served in the q? bins. The covariance matrix due to the Umiss 
fits is determined by 


A S 
C; = | —— ) 3.82 6a PUN)? 10 
J (=) = Eia Eja | ( ebel , ( ) 
where 5(VS...) is the systematic uncertainty of NS, associ- 


ated with the fit to the corresponding Umiss distribution. 


Form factor. To estimate the systematic uncertainty associ- 
ated with the form factor model used to generate signal events 
in the MC simulation, we re-weight the signal MC events so 
that the q? spectra agree with the measured spectra. We then 
re-calculate the partial decay rates in different q? bins with 
the new efficiency matrices which are determined using the 
weighted MC events. The covariance matrix due to this source 
is assigned via C;; = 0(AT';)ó(AT ;), where 6(AT;) denotes 
the change of the partial width measurement in the ith q? bin. 


FSR recovery. To estimate the systematic covariance matrix 
associated with the FSR recovery of the positron momentum, 
we remeasure the partial decay rates without the FSR recov- 
ery. The covariance matrix due to this source is assigned via 
Ci; = Ó6(AT;)Ó(AT ;), where 9(AT';) denotes the change of 
the partial decay rate measurement in the ith q? bin. 


MC statistics. The systematic uncertainties due to the lim- 
ited size of the MC samples used to determine the efficiency 


10 


TABLE VI. Summary of the measured partial decay rates, relative statistical uncertainties, systematic uncertainties and corresponding corre- 


lation matrices for Dt — KP? e^ v, and D* — «?e*v.. 
D* + K?e*y, 
q? bin No. 1 2 3 4 5 6 7 8 9 
AT (ns!) 16.97 15.29 13.57 11.65 9.33 7.06 4.96 2.97 1.01 
stat. uncert. (%) 1.45 1.61 1.75 1.91 2.12 2.44 2.92 3.77 6.56 
stat. correl. 1.000 
—0.073 1.000 
0.001 —0.084 1.000 
0.000 0.003 —0.091 1.000 
0.000 0.000 0.004 —0.085 1.000 
0.000 0.000 0.000 0.004 —0.085 1.000 
0.000 0.000 0.000 0.000 0.004 —0.075 1.000 
0.000 0.000 0.000 0.000 0.000 0.004 —0.069 1.000 
0.000 0.000 0.000 0.000 0.000 0.000 0.003 —0.059 1.000 
syst. uncert. (%) 3.24 3.10 2.95 2.88 3.02 3.05 2.85 2.54 2.93 
syst. correl. 1.000 
0.981 1.000 
0.979 0.976 1.000 
0.979 0.977 0.973 1.000 
0.978 0.976 0.973 0.970 1.000 
0.974 0.972 0.970 0.970 0.965 1.000 
0.966 0.964 0.963 0.962 0.960 0.954 1.000 
0.932 0.930 0.929 0.929 0.926 0.923 0.911 1.000 
0.891 0.889 0.886 0.888 0.886 0.883 0.875 0.840 1.000 
D? — «ety, 
qd? bin No. 1 2 3 4 5 6 7 
AT (ns!) 0.664 0.578 0.474 0.477 0.432 0.503 0.372 
stat. uncert. (%) 4.55 5.53 6.60 6.48 7.28 6.52 8.97 
stat. correl. 1.000 
—0.122 1.000 
0.011 —0.171 1.000 
—0.002 0.019 —0.190 1.000 
0.000 —0.003 0.021 —0.190 1.000 
0.000 —0.001 —0.005 0.016 —0.167 1.000 
—0.002 —0.003 —0.003 —0.008 —0.004 —0.128 1.000 
syst. uncert. (%) 1.53 1.52 1.51 1.61 1.88 1.92 1.73 
Syst. correl. 1.000 
0.739 1.000 
0.742 0.664 1.000 
0.758 0.737 0.650 1.000 
0.772 0.740 0.712 0.698 1.000 
0.781 0.749 0.711 0.760 0.772 1.000 
0.760 0.730 0.697 0.727 0.756 0.740 1.000 
matrices are translated to the covariance via The total systematic covariance matrix is obtained by sum- 
i > ming all these matrices. Table VI summarizes the relative 
ye a yê -1 .-1 size of systematic uncertainties and the corresponding corre- 
Cis (— x] 2 (NS, Nas coles Sg) QD lations in the measurements for the partial decay rates of the 


aß 


where the covariance of the inverse efficiency matrix elements 
are given by [18] 


cov[ezs ea] = J (€ 


ij 


Nor (es? (ee) 2) 


K? and n° decay branching fractions. The systematic un- 
certainties due to the branching fractions of K? > ntn 
(0.07%) and 7° — yy (0.03%) are fully correlated between 

2 . 
q* bins. 


D+ — K?etve and D+ — n?etve semileptonic decays. 


VI. FORM FACTORS 


To determine the product f+ (0)|V.s(a)| and other form fac- 
tor parameters, we fit the measured partial decay rates using 
Eq. (1) with the parameterization of the form factor f, (q?). 
In this analysis, we use several forms of the form factor pa- 
rameterizations which are reviewed in Sec. VI A. 
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A. Form factor parameterizations 


In general, the single pole model is the simplest approach 
to describe the q? dependence of the form factor. The single 
pole model is expressed as 


__ £40) 
fl) Pm P/M i 


(13) 


where f4(0) is the value of the form factor at q? = 0, and 
Mpole i$ the pole mass, which is often treated as a free param- 
eter to improve fit quality. 

The modified pole model [19] is also widely used in Lat- 
tice QCD (LQCD) calculations and experimental studies of 
these decays. In this parameterization, the form factor of the 
semileptonic D — Pe* v, decays is written as 


f+(0) 


SEEE A 14 
HOS cpm a C^ 


where m p»+ is the mass of the D meson, and a is a free 


(5) 
parameter to be fitted. 
The /SGW2 model [20] assumes 


2 —2 
P = n (14 lla) . a» 


where q2 ax is the kinematical limit of g?, and r is the conven- 
tional radius of the meson. 

The most general parameterization of the form factor is the 
series expansion [21], which is based on analyticity and uni- 
tarity. In this parameterization, the variable q? is mapped to a 
new variable z through 


te= ei 
m, (16) 
t4 — g + ytt to 

with t4. = (mp+ + mp)? and to = ty (1 — y1 -— t. /t,). 
The form factor is then expressed in terms of the new variable 
zas 


au. to) ES 


HE M m 
P(@)0(@, to) 2. k(to)[z(a^,t0)]^, aD 


where aj(to) are real coefficients. The function P(g?) is 
P(g?) = z(t, mh) for D > K and P(g?) = 1 for D > m. 
The standard choice of $(q?, to) is 


2 gy — (m8 VT (2t. 0) VP (atto) 7 
«es-() CER) GET) 
4 (sey (t. - d) 
ig (t+ — to)1/4’ 


where me is the mass of the charm quark. 
In practical use, one usually makes a truncation of the above 
series. After optimizing the form factor parameters, we obtain 


f-(0)P(0)9(0, t9) (1 + 32% rx [z(?, to)]^) 
P(o, to (1 + Ye relz(0, to)]*) 


ae = 


(18) 


f 


, 
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where ry = aj (to)/ao(to). In this analysis we fit the mea- 
sured decay rates to the two- or three-parameter series expan- 
sion, i.e., we take kmax = 1 or 2. In fact, the z expansion with 
only a linear term is sufficient to describe the data. Therefore 
we take the two-parameter series expansion as the nominal 


parameterization to determine fE 9o) and |V..,(4)|. 


B. Fitting partial decay rates to extract form factors 


In order to determine the form factor parameters, we fit the 
theoretical parameterizations to the measured partial decay 
rates. Taking into account the correlations of the measured 
partial decay rates among q? bins, the x? to be minimized in 
the fit is defined as 


x? = M (Ar; - AT PCR (AT; — ATH), (20) 


ij 


where AT; is the measured partial decay rate in the ith q? bin, 
C TA is the inverse matrix of the covariance matrix C;;. In the 
ith q? bin, the theoretical expectation of the partial decay rate 
is obtained by integrating Eq. (1), 


2 
Imax, i G2 
art = [77 xara, QD 
a 


min, 


where Q2 ini and VE are the lower and upper boundaries of 
that q? bin, respectively. 

In the fits, all parameters of the form factor parameteri- 
zations are left free. The central values of the form factor 
parameters are taken from the results obtained by fitting the 
data with the combined statistical and systematic covariance 
matrix together. The quadratic difference between the uncer- 
tainties of the fit parameters obtained from the fits with the 
combined covariance matrix and the uncertainties of the fit pa- 
rameters obtained from the fits with the statistical covariance 
matrix only is taken as the systematic error of the measured 
form factor parameter. The results of these fits are summa- 
rized in Table VII, where the first errors are statistical and the 
second systematic. 

Figure 5 shows the fits to the measured differential decay 
rates for Dt + K?e*v, and D+ — s? e*v,. Figure 6 shows 
the projection of the fits onto f4 (q?) for the D+ + K?e*v, 
and Dt — m?e*vw, decays, respectively. In these two fig- 
ures, the dots with error bars show the measured values of the 
form factors, f. (q?), in the center of each q? bin, which are 
obtained with 


AT; 2473 
Die a (22) 
di XGp; Weal 
in which 
2 
Imax i 
2 hu PEP Pa? 
Di ; (23) 


B f+ (as Pes Ginin,i) 
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TABLE VII. Summary of results of form factor fits for D+ + KPetw, and Dt — Tete, where the first errors are statistical and the 


second systematic. 


Decay mode 
Dt + K° ety. 


Dt — Wetv. 


Decay mode 
Dt > K?e*y, 


Dt — metve 


Decay mode 
Dt > K9?e*y, 


Dt — met ve 


Decay mode 
D* 


Decay mode 
D+ > K?ety, 


Dt > a?e*v, 


Single pole model 


f+(0)|Veq 
0.7094 + 0.0035 + 0.0111 
0.1429 + 0.0020 + 0.0009 


f+(0)|Veq 
0.7052 + 0.0038 + 0.0112 
0.1400 + 0.0024 + 0.0010 


f+ (0) |Vea 
0.7039 + 0.0037 + 0.0111 
0.1381 + 0.0023 + 0.0007 


f+ (0) |Vea 
0.7053 + 0.0040 + 0.0112 
0.1400 + 0.0026 + 0.0007 


f+ (0) |Vea 
0.6983 + 0.0056 + 0.0112 
0.1413 + 0.0035 + 0.0012 


p 


IN 
o 
T T T I T T T T T T T [ T T T T 


I ] | 
D*>K’e'v, | 
—e— Data | 
— — Single Pole Model | 
Modified Pole Model4 
---- ISGW2 Model 
Z series (2 par.) 
---- z series (3 par.) 


0.5 1 
a? (GeV?/c^ 


Modified pole model 


ISGW2 model 


Three-parameter series expansion 


1.5 2 


Mpole (GeV/c?) 
1.935 + 0.017 + 0.006 
1.898 + 0.020 + 0.003 


a 
0.294 + 0.031 + 0.010 
0.285 + 0.057 + 0.010 


r (GeV-1c?) 
1.587 + 0.023 + 0.007 
2.078 + 0.067 + 0.011 


Two-parameter series expansion 


Ti 
—2.18 + 0.14 + 0.05 
—2.01 + 0.13 + 0.02 


Ti T2 
—1.76 + 0.25 + 0.06 —13.4 + 6.3 + 1.4 
—2.23 + 0.42 + 0.06 1.4 + 2.5 + 0.4 


2.5 T l 
D'—1?e*v, 
RN —e— Data | 
$ 20 — — Single Pole Model | 
QN |. TN 000€ Modified Pole Model 
> F ---- ISGW2 Model J 
5 15- ] z series (2 par.) 
> E Nye series (3 par.) 
n E Ni " 4 
c 
—10- 
is; L 
= 
s 0.5 - 
0.0 | 
0 1 2 3 


a? (GeV?/c^ 


FIG. 5. Differential decay rates for Dt — K?e*w, (left) and Dt > r°etve (right) as a function of g. The dots with error bars show the 
data and the lines give the best fits to the data with different form factor parameterizations. 


where |V.,| = 0.97351 4 


- 0.00013 and |V,4| = 0.22492 + 


0.00050 are taken from the SM constraint fit [8]. In the cal- 


culation of pj, f (q?) is computed using the two parameter 


series parameterization with the measured parameters. 


C. Determinations of /7 (0) and /7 (0) 


Using the jim (0)|V..(a)| values from the two-parameter 
series expansion fits and taking the values of |Vos(a) | from the 
SM constraint fit [8] as inputs, we obtain the form factors 


K 
+ 


0.004 


0.012 (24) 


(0) = 0.725 
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3.0 : | 
D*—>n etv, yi 
| — Data n | 
2.5 — — — Single Pole Model pan 
spi Modified Pole Model I / 
NEL ISGW2 Model / 4 
ies (2 par: t 
p 20- z series (2 par.) 4 | 
d ---- z series (3 par.) / 
0.5 | | 
0 1 2 3 


q? (GeV?/c^) 


FIG. 6. Projections on f+ (q) for Dt + K° ety, (left) and D+ — n?etve (right) as function of q?, where the dots with error bars show 
the data and the lines give the best fits to the data with different form factor parameterizations. 


and 


f1(0) = 0.622 


0.012 


0.003, (25) 


where the first errors are statistical and the second systematic. 


VII. DETERMINATIONS OF |V..| AND |Vea| 


Using the values of p (0)|Vesca)| from the two- 
parameter z-series expansion fits and in conjunction with the 
form factor values f¥ (0) = 0.747 + 0.011 + 0.015 [22] and 


f£(0) = 0.666 + 0.020 + 0.021 [23] calculated from LQCD, 
we obtain 

|Ves| = 0.944 + 0.005 + 0.015 + 0.024 (26) 
and 

|V.a| = 0.210 + 0.004 + 0.001 + 0.009, (27) 


where the first uncertainties are statistical, the second system- 
atic, and the third are due to the theoretical uncertainties in the 
LQCD calculations of the form factors. 


VIII. SUMMARY 


In summary, by analyzing 2.93 fb^! of data collected 
at 3.773 GeV with the BESIII detector at the BEPCII, the 
semileptonic decays for D^ > K?e*v, and D+ + n?etve 
have been studied. From a total of 1703054 + 3405 D- tags, 
26008 + 168 Dt + Ke*yv, and 3402+ 70 Dt > s? e*v, 
signal events are observed in the system recoiling against the 
D- tags. These yield the absolute decay branching fractions 
to be B(D* > K?e*v,) = (8.60 + 0.06 £0.15) x 107? and 
B(D* — wet vy.) = (3.63 + 0.08 + 0.05) x 1073. 


We also study the relations between the partial decay rates 
and squared 4-momentum transfer q? for these two decays 
and obtain the parameters of different form factor parame- 
terizations. The products of the form factors and the related 
CKM matrix elements extracted from the two-parameter se- 
ries expansion parameterization are selected as our primary 
results. We obtain f, (0)|V.,| = 0.7053 + 0.0040 + 0.0112 
and f/,(0)|V.4| = 0.1400 + 0.0026 + 0.0007. Using the 
global SM fit values for |V.,| and |V.a|, we obtain the form 
factors ff (0) = 0.725 + 0.004 + 0.012 and f7(0) = 
0.622 + 0.012 + 0.003. Furthermore, using the form factors 
predicted by the LQCD calculations, we obtain the CKM ma- 
trix elements |V.,| = 0.944 + 0.005 + 0.015 + 0.024 and 
|V.a| = 0.210 + 0.004 + 0.001 + 0.009, where the third errors 
are dominated by the theoretical uncertainties in the LOCD 
calculations of the form factors. 
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